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Abstract 

We investigate Laplacians on supercritical bond-percolation graphs with different 
boundary conditions at cluster borders. The integrated density of states of the 
Dirichlet Laplacian is found to exhibit a Lifshits tail at the lower spectral edge, 
while that of the Neumann Laplacian shows a van Hove asymptotics, which results 
from the percolating cluster. At the upper spectral edge, the behaviour is reversed. 
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1 Introduction and summary 



Ever since Mark Kac posed the question "Can one hear the shape of a drum?" 
[15], there has been a great deal of interest in finding relations between the 
geometry of a manifold or a graph and spectral properties of the Laplacian 
defined on it. The impressive works [8-10,6,3], which have been chosen by way 
of example, witness the steady progress achieved in recent years and provide 
further references. Whereas Laplacians on manifolds dominated the scene in 
the earlier years, the rise of spectral graph theory [23,22,11,5,7] in the late 
1980s and 90s has contributed to deepen our understanding of the discrete 
case. 
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Spectral theory of random graphs, however, is still a widely open field. The 
very recent contributions [21,2,14] take a probabilistic point of view to de- 
rive heat-kernel estimates for Laplacians on supercritical Bernoulli bond- 
percolation graphs in the rf-dimensional hyper-cubic lattice. On the other 
hand, traditional methods from spectral theory are used in [19] to investigate 
the integrated density of states of Laplacians on subcritical bond-percolation 
graphs. Depending on the boundary condition that is chosen at cluster bor- 
ders, two different types of Lifshits asymptotics at spectral edges were found 
[19]. For example, the integrated density of states of the Neumann Laplacian 
behaves as 

" N N (E) - iV N (0) ~ exp{-£- 1/2 } " as E j (1.1) 

at the lower spectral edge for bond probabilities p below the percolation thresh- 
old p c . We have put quotation marks here, because, strictly speaking, one 
should take appropriate logarithms on both sides. The Lifshits exponent 1/2 
in (1.1) is independent of the spatial dimension d. This was explained by the 
fact that, asymptotically, N-^ is dominated by the smallest eigenvalues which 
arise from very long linear clusters in this case. In contrast, for the Dirichlet 
Laplacian and p < p c , it was found that 

" N B {E) ~ exp{-E' d/2 } " as E 10. (1.2) 

We note that iV D (0) = 0. The Lifshits exponent in (1.2) comes out as d/2, 
because the dominating small Dirichlet eigenvalues arise from large fully con- 
nected cube- or sphere-like clusters. Thus, depending on the boundary condi- 
tion (and the spectral edge) different geometric graph properties show up in 
the integrated density of states. We refer to the literature cited in [19] for a 
discussion of other spectral properties of these and closely related operators, 
for the history of the problem and what is known in the physics literature. 
Lifshits asymptotics for a Neumann Laplacian on Erdos-Renyi random graphs 
are studied in [18]. 

In this paper we pursue the investigations of [19] and ask what happens to 
(1.1) and (1.2) in the supercritical phase of bond-percolation graphs. Clearly, 
one would not expect the contribution of the finite clusters to alter the picture 
completely. But for the infinite percolating cluster, the story may be differ- 
ent. Indeed, we will prove that the percolating cluster produces a van Hove 
asymptotics 

" N N (E) - JV N (0) ~ E d/2 " as E | (1.3) 

in the Neumann case for p > p c . There is also an additional Lifshits-tail 
behaviour due to finite clusters, but it is hidden under the dominating asymp- 
totics (1.3). Loosely speaking, (1.3) is true because the percolating cluster 
looks like the full regular lattice on very large length scales (bigger than the 
correlation length) for p > p c . On smaller scales its structure is more like that 
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of a jagged fractal. The Neumann Laplacian does not care about these small- 
scale holes, however. All that is needed for (1.3) to be true is the existence 
of a suitable rf-dimensional, infinite grid. In contrast, the Dirichlet Laplacian 
does care about holes at all scales so that (1.2) continues to hold for p > p C) as 
we shall prove. Low-lying Dirichlet eigenvalues require large fully connected 
cube- or sphere-like regions, and this is a large-deviation event. 

Closely related large-deviation results for Laplacians on percolation graphs 
have been obtained in [1,4]. To be precise, [1,4] refer to the Pseudo-Dirichlet 
Laplacian Ag in the sense of our Definition 2.1(h) below. Considering both 
site- and bond-percolation graphs, and using a discrete version of the method 
of enlargement of obstacles, Antal [1] derives the long-time asymptotics for 
the mean (i.e. annealed) hitting-time distribution of the set of absent sites 
(resp. bonds) for the random walk generated by Ag. Biskup and Konig work 
in the setting of the parabolic Anderson model, which contains Ag on site- 
percolation graphs as a special case. In particular, they establish a Lifshits tail 
for the corresponding integrated density of states, see also Remark 2.6 (v). 

This paper is organised as follows. In the next section we give a precise state- 
ment of our results in Theorems 2.5 and 2.7. Section 3 is devoted to the proof 
of Theorem 2.5. In this proof we follow the strategy laid down in [24], see also 
[25]. The goal there was to establish Lifshits tails in the context of random 
Schrodinger operators. Finally, Section 4 contains the proof of Theorem 2.7, 
where we apply the recent deep heat-kernel estimates from [21,2,14]. 



2 Definitions and precise formulations 



To set up the mathematical arena, let us first recall some notions from 
Bernoulli bond percolation. For d e N , a natural number, we denote by h d 
the (simple hyper-cubic) lattice in d dimensions. Being a graph, the lattice 
h d = (Z d , E d ) has the vertex set Z rf and the edge set K d given by all unordered 
pairs {x, y} of nearest-neighbour vertices x, y G Z d , that is, those vertices 

/ d „\ 1/2 

which have Euclidean distance \x — y\ := [J2 u =i \x v — y u \ ) =1- Here, el- 
ements of Z d are canonically represented as ci-tuples x = (xi, . . . , xj) with 
entries from Z. Next, we consider the probability space = {0, l} Ed , which 
is endowed with the usual product sigma-algebra, generated by finite cylinder 
sets, and equipped with a product probability measure P. Elementary events 
in Q are sequences of the form u = (uj{ Xty y){ x ^ e ^d, and we assume their en- 
tries to be independently and identically distributed according to a Bernoulli 
law P(u{ x ^} = 1) = p with bond probability p e]0, 1[. To a given u e £1, we 
associate an edge set := <{x, y} G E d : 0J{ x , y } — l}- 
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A bond-percolation graph in Z d is the mapping W : f2 3 u i— > := (Z d , 
with values in the set of subgraphs of h d . Given x G Z d , the vertex degree 
dcg(u) (x) counts the number of edges in which share x as a common vertex. 



Definition 2.1. The random operators D : VI 3 uo i— > D^, respectively 
A : 11 9 w h are called vertex-degree operator, respectively adja- 

cency operator, of bond-percolation graphs in 7L d . Their realisations, : 
£ 2 (Z d ) -> £ 2 (Z d ), respectively : £ 2 (Z d ) -> £ 2 (Z d ), act on the Hilbert space 
of complex-valued, square-summable sequences indexed by Z d according to 

D^(p{x) : = dcg(w) (x) ip(x) , 
?/ez d : {x,»}e<?( w ) 

for all G £ 2 (Z d ), all x e Z d and all w G 11. With these definitions, we 
introduce Laplacians on bond-percolation graphs for three different "boundary 
conditions" at non-fully connected vertices 

(i) Neumann Laplacian: An := D — A, 

(ii) Pseudo-Dirichlet Laplacian: Ag := An + (2dll — D) = 2d\l — A, 

(iii) Dirichlet Laplacian: A D := A N + 2(2dll — D). 

Here 11 stands for the identity operator on £ 2 (Z d ). 

Remarks 2.2. (i) The motivation and origin of the terminology for 

the different boundary conditions are discussed in [19] - together with some 
basic properties of the operators. 

(ii) The random self-adjoint Laplacians are ergodic with respect to 
Z d -translations. Hence, their spectra and the spectral subsets arising in the 
Lebesgue decomposition are all equal to non-random sets with probability one. 
In particular, the spectrum is P-almost surely given by spec(Ax) = [0, Ad] for 
all X G {N, D, D}, as was shown in [19]. 

Next, we define the quantity of our main interest for this paper, the integrated 
density of states of A x . To this end let 5 X G £ 2 (Z d ) be the sequence which is 
concentrated at the point x G Z d , i.e. S x (x) := 1 and 5 x (y) := for all 
y G Z d \ {x}. Moreover, O stands for the Heaviside unit-step function, which 
we choose to be right continuous, viz. Q(E) := for all real E < and 
Q(E) := 1 for all real E ^ 0. 

Definition 2.3. For every p G]0,1[ and every X G {N,D,D} we call the 
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function 



N x : R 3 E ^ N X {E) := ^P(dcu) {6 ,G(E - A^)S ) (2.2) 

with values in the interval [0, 1] the integrated density of states of A x . 

Remarks 2.4. (i) The integrated density of states Nx is the right- 

continuous distribution function of a probability measure on R. The set of its 
growth points coincides with the P-almost-sure spectrum [0, Ad] of Ax- 

(ii) It is shown in [19] that the Laplacians are related to each other by 
a unitary involution, which implies the symmetries 

iV D(N) (£) = l- lim N m) (e) 

£ 1 4 a — a, 

for their integrated densities of states for all The limits on the right- 

hand sides of (2.3) ensure that the discontinuity points of Nx are approached 
from the correct side. 

(iii) By ergodicity, Definition 2.3 of the integrated density of states co- 
incides with the usual one in terms of a macroscopic limit of a finite-volume 
eigenvalue counting function. More precisely, let A C Z d stand for bounded 
cubes centred at the origin with volume |A|. For every X e {N, D, D} let Ax,a 
be the finite- volume restriction of A x to £ 2 (A) introduced in Def. 1.11 in [19]. 
Then there exists a set il' C f2 of full probability, P(f2') = 1, such that 



Nx(E) = lim 



^ trace^e^-A^) 



(2.4) 



holds for all wGfl' and all E e R , except for the (at most countably many) 
discontinuity points of Nx, see Lemma 1.12 in [19]. In Section 3 we will con- 
struct another finite-volume restriction of Ag, for which (2.4) holds, too. 

Let p c = p c (d) denote the critical bond probability of the percolation transition 
in Z d . We recall that p c — 1 for d — 1, otherwise p c e]0, 1[, see e.g. [13]. Despite 
the title of this paper, our first main result covers the non-percolating phase 
p E]0,p c [ and the critical point p = p c , too. 

Theorem 2.5. Assume d e N and p e]0, 1[. Then the integrated density of 
states Nx of the Laplacian A x on bond-percolation graphs in Z d exhibits a 
Lifshits tail at the lower spectral edge 

Hhyi = _g for Xe{B , D} (2 . 5) 

eio InE 2 J 1 s y ' 
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and at the upper spectral edge 



^ h, ^: N i m =-j /- X e{ N,D}. (2.6) 

Remarks 2.6. (i) The theorem follows directly from the upper and 

lower bounds in Lemma 3.1 below, together with the subsequent Re- 
mark 2.6 (ii). In fact, the bounds of Lemma 3.1 provide a slightly stronger 
statement than Theorem 2.5. 

(ii) The Lifshits tails at the upper spectral edge are related to the ones 
at the lower spectral edge by the symmetries (2.3). 

(iii) In the non-percolating phase, p e]0,p c [, the content of the theorem 
is known from [19], where it is proved by a different method. The method of 
[19], however, does not seem to extend to the critical point or the percolating 
phase, p e]p c , 1[. 

(iv) The Lifshits asymptotics of Theorem 2.5 are determined by those 
parts of the percolation graphs, which contain large, fully-connected cubes. 
This also explains why the spatial dimension enters the Lifshits exponent d/2. 

(v) We expect that (2.5) can be refined in the case X = D as to obtain 
the constant 

In N?.(E) 

1- -^==-*(M- (2-7) 

An analogous statement is known from Thm. 1.3 in [4] for the case of site- 
percolation graphs. Moreover, it is demonstrated in [1] that the bond- and the 
site-percolation cases have similar large-deviation properties. 

Our second main result complements Theorem 2.5 in the percolating phase. 

Theorem 2.7. Assume d£N \ {1} and p e]p c , 1[. Then the integrated den- 
sity of states of the Neumann Laplacian A N on bond-percolation graphs in Z d 
exhibits a van Hove asymptotics at the lower spectral edge 

eio \nE 2 y J 

while that of the Dirichlet Laplacian Ad exhibits one at the upper spectral edge 

ln{N»(4d)-N D (E)] _d 

ln(4d - E) " 2 ' (2 ' 9) 

where N D {4d) := lim EUd N D {E) = 1 - JV N (0). 

Remarks 2.8. (i) The theorem follows directly from the upper and 

lower bounds in Lemma 4.1 below, together with the symmetries (2.3). In 
fact, the bounds of Lemma 4.1 provide a slightly stronger statement than 
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Theorem 2.7. Lemma 4.1 relies mainly on recent estimates [21,2,14] for the 
long-time decay of the heat kernel of An on the percolating cluster. 

(ii) The reference value A^n(O) in (2.8) results from the mean number 
density of zero eigenvalues of the Neumann Laplacian [19]. It is given by 

N N (0) = K(p) + (l-p) 2d , (2.10) 

where n(p) is the mean number density of clusters, see e.g. Chap. 4 in [13], 
and (1 - p) 2d the mean number density of isolated vertices. 

(iii) The counterpart of Theorem 2.7 for the non-percolating phase, 
p e]0,p c [, was proved in [19]. There, iV N was shown to have a different kind of 
Lifshits asymptotics with a Lifshits exponent 1/2 at the lower spectral edge, 
see also Section 1, and the same is true for Nn at the upper spectral edge. 
This type of Lifshits behaviour is caused by large (isolated) linear clusters, 
which explains why the spatial dimension does not influence the Lifshits ex- 
ponent. This behaviour is also present for p e]p c , 1[, but hidden under the 
more dominant van Hove asymptotics caused by the percolating cluster. 

(iv) At the critical point p = p c , the behaviour of N-^ at the lower 
spectral edge, respectively that of iVb at the upper spectral edge, is an open 
problem. 



3 Proof of Theorem 2.5 

In this section we prove the Lifshits-tail behaviour of Theorem 2.5. Thanks to 
the symmetries (2.3), it suffices to consider the lower spectral edge only. 

Lemma 3.1. For every d G N and every p G]0,1[ there exist constants er>, 
a u , ai g]0,oo[ such that 

exp{- ai E~ d / 2 } ^ N D (E) ^ N^(E) <: e^{-a u E~ d / 2 } (3.1) 

holds for all E g]0,£d[- 

Proof. The left inequality in (3.1), i.e. the lower bound on iV D , was proved 
in Lemma 2.9 in [19]. The middle one simply reflects the operator inequality 
^ \ which is valid for all u G Vt. So it remains to prove the upper 
bound on N~. 

We follow the strategy of the proof in [24], see also Chap. 2.1 in [25]. To do so, 
we have to fix some notation, first. Given a bounded cube A C Z d and x G A, 
we introduce the boundary degree 

b dA (x):=\{{x,y}eE d :y<£A}\ (3.2) 
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as the cardinality of the set of edges in the regular lattice h d that connect x 
with Z d \ A. The restriction := (A, S^) with := {{x,y} G <S^ : 

x, y G A| of any realisation of a bond-percolation graph to A is obtained 
by keeping only vertices and edges within A, and cL( w ) (x) ^ 2d — bdA(x) stands 

for the associated vertex degree of x G A. In particular, E A := G E d : 

x,y £ A| is the edge set of the fully connected cube L A := (A, E A ), that is the 

restriction of the regular lattice h d to A. Finally, let £ 2 (A) be the Hilbert space 
of complex-valued (finite) sequences indexed by A, and, given any subgraph 
<& := (A, (£) of L A , we introduce the operator : £ 2 (A) — > I 2 (A), (p i— > fj©y?, 
where 

:= - E ^(y) + (2d-b dA (x)^ip(x) 

yeA:{x,y}e<£ 

= ]T {<P(x) - <p(y)) + (2d - b 9A (x) - d e (x))<p(x) (3.3) 

ye\:{x,y}e<£ 

for all if G £ 2 (A) and all x G A. Now, we define the restriction of the Pseudo- 
Dirichlet Laplacian Ag to the cube A with Neumann conditions along the 
boundary <9A of A as the random bounded self- adjoint operator H A with re- 
alisations := for all uo G Vt. 

Next we claim that 

(3.4) 



N$(E) = inf 

D AcZ d 



i- jf P(da/) trace, 2(A) G(E - H^) 



holds for all £eR. This is so, because (i) the operator differs from 
the finite-volume restriction A~^ in Remark 2.4 (iii) by a perturbation whose 
rank is at most of the order of \dA\, the surface area of the cube A. Hence, 
(2.4) remains true for X = D and with A~^ being replaced by on its 

right-hand side, 



N%(E) = lim 



i- trace, 2(A) Q(E - H^) 



(3.5) 



On the other hand, H^ff* is designed in such a way that H$ © 



h a%a 2 nolds on ^ 2 (AiUA 2 ) for all bounded cubes A : , A 2 C Z d with AiHA 2 = 
and for all oo G ft. Hence, Q{E-H { A ] ) gives rise to a subergodic process and we 
conclude from the Ackoglu-Krengel subergodic theorem that the right-hand 
side of (3.5) equals the right-hand side of (3.4) - again for all continuity points 
of the limit and uniformly for u in a set of probability one. (iii) From this 
we have (3.4) for all continuity points of both sides. But since both sides of 
(3.4) are right-continuous functions of E, equality holds for all E G K, and 
the derivation is complete. 



8 



From (3.4) we infer the upper bound 



N^(E)<: inf F[E A ^ E] , (3.6) 

where the non-negative random variable E A stands for the smallest eigenvalue 
of the random operator H\. 

The aim is to obtain a simple large-deviation estimate for the probability 
in (3.6). This will be achieved with the help of analytic perturbation theory 
along the lines of [24], see also Sec. 4.1.10 in [25]. We write H 0A := $j h d for 
the Neumann Laplacian of the fully connected cube L A and W A ■= Ha — H 0:A . 
Given t G [0, 1], we introduce 

H A (t) :=H 0A + tW A (3.7) 

so that H A (1) = Ha- We want to construct an upper bound for the probability 
that Ea is small. Denoting the bottom eigenvalue of H A {t) by E A (t), we use 
the following ideas. 



(a) The function [0,1] 3 t i— > E A (t) is non-decreasing, -E'a(O) = and 
E A {l)=E A . 

(b) This function can be linearised, if its argument is small enough. More 
precisely, there exist constants r, f3 e]0, oo[, which depend only on the 
spatial dimension d, such that 

\E A (t)-tE' A (0)\^(3t 2 \A\ 2 / d (3.8) 

for all t G [0,r\A\- 2 / d ]. Here, we have set E' A {0) := ^.E7 A (t)| 

(c) The slope E' A (0) obeys a large-deviation estimate. Given any a e]0, 1— p[, 
there exists a constant 7 e]0, oof, which depends on p and d, such that 

F[E' A (0) <:a]^ e^ |A| . (3.9) 



We will prove (a) with a Perron-Frobenius argument in Lemma 3.2 and discuss 
observations (b) and (c) below. For the time being, let us go on to estimate 
the probability that E A is small. 

Suppose E A (t) ^ E. Then we conclude from (a), the triangle inequality and 
(b) that 



E' A (0) < ^# + 



E A (t) 



- E' A (0) 



E 



+ (3t\K\ 2/d , 



(3.10) 



provided t is small enough. So we need to adjust t = and A = such 
that tE ^ r\A E \~ 2 / d . Moreover, we aim to achieve that the right-hand side of 
(3.10) is bounded from above by some a < 1 — p. This is accomplished in the 
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following way. Without restriction we can assume that, in addition, a < 2(3t. 
Then we set t E ■— a/(2(3\A E \ 2 ^ d ) and choose the size of the cube such that 



a 



- 1 < |A £ | lAi < 



a 



(3.11) 



2( / 3£) 1 /2 



2((3Ey/ 2 ' 



For this to make sense, the right-hand side of (3.11) has to exceed 2. So, we 
restrict ourselves to low energies, say E e]0, £d[, and summarise this argument 

as 



Note that e-q depends only on p and d. 

Altogether, we infer from Eq. (3.6), observation (a), implication (3.12) and 
observation (c) that 



N^E) ^ ¥[E AE (t E ) ^E]^ ¥[E' Ae (0) < a] < e'^ ^ e -^~ d/2 , ( 3 .13) 



where a u e]0, oo[ is a constant that depends only on p and d. 

Next, we verify observations (b) and (c) above. Observation (b) relies on a 
deterministic result from analytic perturbation theory. To this end we consider 
the operator family H(z) := H + zEL\ for z G C. Here, H := H 0i a is the 
Neumann Laplacian of and Hi := the perturbation with w G O 

arbitrary, but fixed. The bottom eigenvalue of H is an isolated simple 
eigenvalue. Its isolation distance "d := dist^O, spec(i/o\ {0})) is determined by 
the magnitude of the smallest non-zero eigenvalue of H . This distance obeys 
the estimate d ^ c\A\~ 2 ^ d for some constant c e]0,oo[, which follows from 
reducing the eigenvalue problem for H to that of a linear chain by separation 
of variables and applying a Cheeger-type inequality, see e.g. (2.6) in [19]. 
Moreover, we have the uniform bound ||W^|| ^ 2d for the operator norm of 
the perturbation so that H(z) has one isolated eigenvalue E(z) in the complex 
disc B# /2 (0) provided z < #/(4d). We refer to [17], Chap. II, §1, Sees. 1, 2 
and Chap. VII, §3, Sees. 1, 2 and 4 for a detailed exposition of the general 
method. Elementary function theory then gives an estimate for the second 
derivative of E(z), and Taylor's theorem yields (b). Details of the argument, 
geared towards our application here, can also be found in Sec. 4.1.10 in [25]. 

Concerning observation (c), we refer again to analytic perturbation theory. 
The Feynman-Hellmann formula yields 



where </?o := |A| 1 ^ 2 , the normalised vector in £ 2 (A) with equal components, 
is the ground state of the unperturbed operator H A . Therefore, recalling 



E Ae (t E ) ^ E implies E' Ae (0) ^ a < 1 - p . 



(3.12) 



E' A (0) = (<p ,W A <po) , 



(3.14) 
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Wa = Ha — if ,A = $)y A — an d the definition in (3.3), we have 

4 w), (o) = ^ E E i = 4, E 

^ m^T E ( x -"{*,!/}) ( 3 - 15 ) 

™ {x,y}eEf 

for all w G fi. We recall that the w/^i's, which indicate the presence of 
an edge in the bond-percolation graph, are i.i.d. Bernoulli distributed with 
mean p. Hence, (3.9) follows from standard large-deviation estimates, see e.g. 
inequality (27.4) in [16] or Thm. 1.4 in [26]. □ 

So far we have deferred the proof of observation (a) in the above demonstra- 
tion. This is a deterministic result which we address now in 

Lemma 3.2. Let A C Z d be a bounded cube, let (25 = (A, (£) be a subgraph of 
the fully connected cube and let ft® be the finite-volume Laplacian (3.3) on 
£ 2 (A). Fort G R let t(t) be the smallest eigenvalue of 

f)(t) :=S) hi +tW, (3.16) 
where 2U := i}® — f) L d . Then the function [0, 1] 3 1 1— > e(i) is non- decreasing. 

Proof. We observe from the definition of 2U and (3.3) that 

2Uy?(x) = E (3.17) 

j/€A:{s,j/}6E£\<£ 

for all if G £ 2 {A) and all x G A. Given i G [0, 1] let us rewrite t)(t) = - (1 - 
t) 2H =: 2dll — a(t). In particular, o(l) = 2dll — S)® acts as 

a(l)(p(x)= E <p(y) +b d A(x)<p(x) (3.18) 

S/eA:{a:,s/}ee 

for all if G £ 2 (A) and all x G A. Equations (3.17) and (3.18) show that the 
self- adjoint linear operator 

a{t) = a(l) + (l-t)2U, (3.19) 

which lives on the finite-dimensional Hilbert space £ 2 (A), has only non- 
negative matrix elements (8 X , a(t)S y ) for all x,y G A. Together with the min- 
max principle, this implies that one can choose the eigenvector (s) correspond- 
ing to the largest eigenvalue of a(t) in such a way that all their components 
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in the basis {S x } x( za are non- negative. Hence, the same is true for the eigen- 
vector (s) corresponding to the smallest eigenvalue of f)(t). Thus, another ap- 
plication of the min-max principle yields 

e (t 2 )= inf MjhM - inf M(tM + & - h^Wy) 

0^e£ 2 (A) (ip,tp) O^G^(A) (<p,(p) 

> t(h) (3.20) 

for all ^ ti ^ t 2 ^ 1, because the scalar product involving 2U is non- negative 
by (3.17). □ 



4 Proof of Theorem 2.7 



In this section we prove the van Hove asymptotics of Theorem 2.7. Again, it 
suffices to consider the lower spectral edge, because of the symmetries (2.3). 
That asymptotics follows from 

Lemma 4.1. Assume d G N \ {1} and p G]p c , 1[. Then there exist constants 
£n> C u , Ci g]0, oo[ such that 

QE d / 2 ^ N N (E) - JV N (0) < C u E d ' 2 (4.1) 

holds for all E g]0, £n[- 

To prove Lemma 4.1 we separate the contribution of the percolating cluster 
to iV N from that of the finite clusters. 

Definition 4.2. Let f2oo denote the event that the origin belongs to the per- 
colating cluster and, for define 

AW£) := / P(<M (8 ,e(E - aP)6 ) , (4.2) 

which is the contribution of the percolating cluster to the integrated density 
of states of the Neumann Laplacian. We write N Nj00 (t) := J£°dN Nj00 (E) e~ Et 
for its Laplace transform, where t G [0, oo[. 

As is well known, the Laplace transform of (4.2) can be related to the mean 
return probability of a continuous-time, simple random walk {Z t } te [ 0jOO [ on the 
percolating cluster. More precisely, this random walk is the Markov process 
on Z d defined by the following set of rules: Suppose the process is at x G 
Z d . Then, after having waited there for an exponential time of parameter 
one, one of the 2c? neighbours of x in Z d , say y, is chosen at random with 
probability l/(2d). If = 1, then the process jumps immediately to y, 
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otherwise there will be no move. The procedure then starts afresh. Assuming 
that Zq = xq G 1, d is the starting point of the process, we denote its law by 
■ The process Z t is generated by the Neumann Laplacian in the sense 
that the transition probability for going from x to y within time t is given by 

PjtfiZs+t = y\Z s = x) = ^\Z t = y) = (S y , e~ tA ^/^5 x ) for all s G [0, oof 
and all Xo in the same connected component as x and y. From this it follows 
that 

N Nj00 (t) = [ P(dcu) ^ ] {Z 2dt = 0) . (4.3) 

Hence, (P(f2oo)) iV Ni00 (t) is the (conditional) mean return probability at time 
2dt for the process on the percolating cluster for p e]p c , 1]. 

Averaged transition probabilities of Z t or related random walks have recently 
been studied in [21,2,14] with elaborate methods. We state a special case of 
the results as 

Proposition 4.3. Assume deN \ {1} and p E]p c , 1[. Then there exist con- 
stants Q, c u g]0, oof and t g]1, oo[ ; all of which depend only on p and d, such 
that 

Cl r d i 2 ^ jv N)00 (t) < c u r d ' 2 (4.4) 

holds for all t G [t , oo[. 

Remark 4.4. In view of (4.3), the lower bound in Proposition 4.3 can be found 
as Eq. (30) in Appendix D in [21]. That paper also contains a "quenched" 
upper bound, i.e. an upper bound for ^^\Z t = y), which is valid for P-almost 
every oo G JIqo- Unfortunately, it is not clear how to take the probabilistic 
expectation thereof, that is, to get an "annealed" upper bound. On the other 
hand, the authors of [14] prove an annealed upper bound in their Thm. 8.1. But 
this bound includes an additional logarithmic factor. The strongest results, 
both annealed and quenched, are those of Barlow [2], and (4.4) is a special 
case thereof. However, his results apply to a random walk which is generated 
by .D _1 An instead of An. Hence, some additional comments are needed to 
adapt his results, and we address this issue now. 

Proof of Proposition 4.3. According to the preceding remark the proposition 
is established, if we show that Barlow's quenched upper bound for the return 
probability, i.e. the special case x — y of the upper bound in Thm. 2 in [2], 
applies also to the random walk generated by An- Eventually, the upper bound 
in Thm. 2 in [2] is reduced to Prop. 3.1 in [2] via Thm. 1, Prop. 6.1, Thm. 5.7 
and Thm. 3.8 - the latter being nothing but the off-diagonal generalisation of 
Prop. 3.1, so we do not need it here. The reduction does not make use of any 
specific properties of the random walk's generator. Hence, all that remains to 
check is Prop. 3.1. in [2]. It turns out that some of the constants in the proof 
of Prop. 3.1 must be modified for our purpose, but this does not have any 
consequences. In addition, the proof of Prop. 3.1 also requires estimates (1.7) 
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and (1.8) of Lemma 1.1 in [2]. Estimate (1.8) follows from estimate (1.5), as 
is argued in the proof of Lemma 1.1(b) in [2]. This argument applies in our 
situation, too. So, in the end, we have to verify the validity of (1.7) and (1.5) 
in [2] for the random walk generated by A N . 

Estimate (1.5) is seen to hold as a special case of Cor. 11 in [12]. The upper 
bound in estimate (1.7) can be inferred, for example, from Thm. II. 5 in [8], 
which is a general result for ultracontractive Markov semigroups. (Actually, 
we refer to the first of the two theorems with the same number II. 5 in [8].) To 
verify statement (ii) of that theorem, one may use the application following it 
together with the reasoning in Cor. V.2 and its proof with the choice ip(x) = 
x. This choice corresponds to a weak isoperimetric inequality which merely 
reflects that the percolating cluster contains an infinitely long path. From this 
point of view, the weak £ -1 / 2 -decay of the upper bound in estimate (1.7) does 
not come as a surprise. Finally, the lower bound in estimate (1.7) arises solely 
from the uniform growth condition \B(x,R)\ ^ const. R d for the volume of a 
ball around x G 7L d with radius R. Such type of results are well known for 
heat kernels on manifolds and also for discrete-time random walks [20,9,10]. 
To employ them here, we decompose 

0> x {Z t = x) = f; (5 X , K n 5 x ) & x ( there are " ^tempted \ (45) 
^ y jumps up to time t J 

using the stochastic independence of all building blocks of Z t . Here, the con- 
traction K : = 11 — An/ (2d) is the transition matrix of a discrete-time random 
walk on Z d , which controls the directions of the jumps of Z t . The number of 
attempted jumps up to time t is governed by a Poisson distribution with mean 
t. Using this, the lower bound in (1.7) follows from Thm. 3 (ii) in [20] applied 
to K. □ 



In order to apply Proposition 4.3 in the proof of Lemma 4.1, we will use two 
elementary Tauberian inequalities. 

Lemma 4.5. Let n be a positive Borel measure on IR + = [0, oof. Suppose 
there are constants to, 8, q, c u g]0, oof such that the Laplace transform Jl(t) : = 
J R+ n(dE) e~ Et exists for all t e [to, °°[ an d obeys 

cit- 5 < p,(t) ^ c u r s . (4.6) 

Then there exist constants Ci,C u e]0, oo[ such that 

QE* ^n([0,E\) ^C U E S (4.7) 

holds for all E ejO,^ 1 ]- 
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Proof. We express //([O, E\) = J R+ /i(dA) 0(1 — \/E) in terms of the right- 
continuous Heaviside unit-step function and observe the elementary inequality 
e -rx _ e -(r-i) e -x ^ _ x ^ ^ e 1_a: , which is valid for all x G 1R + and all 
t G [l,oo [. The upper bound in (4.7) is obvious now. For the lower bound, 
one has to choose r large enough such that the constant, which arises from 
the application of both estimates in (4.6), is strictly positive. □ 

Proof of Lemma 4.1. We set A^fin : = — ^n,oo for the contribution of the 
finite clusters to the integrated density of states and observe 

N N (E) - N N (0) = [N NM (E) - N NM (0)] + N Nj00 (E) (4.8) 

for all E eR. 

Proposition 4.3 and Lemma 4.5 establish the desired van Hove bounds for 
iV N)00 . Therefore, it remains to show that the finite clusters do not spoil this 
behaviour. Since N N ^ n (E) — AT N)fin (0) ^ for all E G [0, oof, only an appropri- 
ate upper bound is required for the contribution of the finite clusters. We shall 
show in Eq. (4.10) below that Nx,& n (E) — iV Ni fi n (0) obeys even a Lifshits-type 
upper bound, which will then complete the proof. 

The Lifshits behaviour for the contribution of the finite clusters in the perco- 
lating phase arises from the cluster-size distribution - in the same way as it 
was shown to arise in the non-percolating phase in [19]. Indeed, we have 

N NM (E) - 7V N ,fin(0) ^ f{lo G fi fin : |^ H | > (dE)' 1 / 2 } (4.9) 

for all E g]0,oo[. Here f2 fin :— ft\ is the event that the origin belongs 
to a finite cluster, say ^o W \ and I^q^I denotes the number of its vertices. 
Inequality (4.9) follows from repeating the steps that lead to the first inequality 
in Eq. (2.24) in [19] with iV N! fi n instead of iV N . For p > p c the cluster-size 
distribution on the right-hand side of (4.9) decays sub-exponentially according 
to Thm. 8.6.5 in [13] so that we obtain 

N NM (E) - N NM (0) ^ ci exp{-^- { ^ 1)/M } (4.10) 

for all E G]0,oo[ with some constants ci,£ G]0,oo[, which depend only on p 
and d. This completes the proof. □ 
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